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Abstract 

We review some recent results surrounding a general mechanism for producing 
chaotic behavior in periodically-kicked oscillators. The key geometric ideas are 
illustrated via a simple linear shear model. 

Introduction 

This paper reviews some recent results on a topic with a considerable history: the 
periodic forcing of limit cycles. Some 80 years ago, van der Pol and van der Mark 
observed that irregularities developed when certain electrical circuits exhibiting stable 
oscillations were periodically forced [26] . Their work stimulated a number of analytical 
studies; see e.g. O [131 1121 H]. Another classical example of driven oscillators is the 
FitzHugh-Nagumo neuron model [6]; the response of this and other models of biological 
rhythms to external perturbations have been extensively studied (see e.g. [3B]). As 
a topic of mathematical study, the dynamics of forced oscillations is well motivated: 
Oscillatory behavior are ubiquitous in physical, biological, and engineered systems, and 
external forcing, whether artificially applied or as a way to model forces not intrinsic 
to the system, is also commonplace. 

In this article, we are not concerned with modeling specific physical phenomena. 
Instead, we consider a generic dynamical system with a limit cycle, and seek to un- 
derstand its qualitative behavior when the system is periodically disturbed. To limit 
the scope of the problem, we restrict ourselves to periodic kicks, or forcings that are 
turned on for only short durations, leaving the limit cycle ample time to restore itself 
during the relaxation period. We are interested in large-time behavior, particularly in 
questions of stability and chaos. 
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As we will show, one of the properties of the limit cycle that plays a key role in 
determining whether the kicked system is stable or chaotic is shear, by which we refer to 
the differential in speed (or angular velocity) for orbits near the limit cycle. A central 
theme of this article is that under suitable conditions, the impact of a kick can be 
substantially magnified by the underlying shear in the unforced system, leading to the 
formation of horseshoes and strange attractors. This does not always happen, however: 
some limit cycles are more vulnerable, and some types of kicks are more effective than 
others. These ideas are discussed in ED] and [H], the material which forms the 
basis of the present review. 

Even though we seek to provide insight into dynamical mechanisms that operate 
under general conditions, we have found the ideas to be most transparent in a very sim- 
ple linear shear model, to which we will devote a nontrivial part of the paper. Sect. 1 
introduces this example and familiarizes the reader with the various parameters (in- 
cluding the one which measures shear); it also reports on results of a numerical study on 
Lyapunov exponents. Sects. 2 and 3 are organized around explaining these simulation 
results. Along the way, we take the opportunity to review a number of mathematical 
ideas which clearly go beyond this one example. Some of the rigorous results reviewed, 
notably those on SRB measures for a relevant class of strange attractors [281 EB 122] , 
are recent developments. With the main ingredients of the linear shear model and the 
relevant mathematical background in hand, we return to a discussion of general limit 
cycles in the final section. 



1 Increasing Shear as a Route to Chaos 

This section introduces the main example we use in this review, and acquaints the 
reader with the various parameters in the model and how they impact the dynamics. 
Of particular interest to us is the effect of increasing shear. Numerically computed 
Lyapunov exponents as functions of shear are presented in Sect. 1.3. They will serve 
as a focal point for some discussions to follow. 

1.1 Periodic kicking of linear shear flow 

Our main example is the periodic kicking of a linear shear flow with a hyperbolic 
limit cycle. This 2D model was studied rigorously in [291 EO] and numerically in [HI 
|35] . Though exceedingly simple in appearance, it already exhibits rich and complex 
dynamical behaviors. 
The model is given by 

e = 1 + cry , 

y = -\y + A-smi2ne)-Zn=-oo^it-^^) ^ 

where {9, y) eS^ x M, S-*^ = R/Z, and a, A, A and r are constants with a. A, r > 0. We 
will refer to Eq. ([T]) with A = as the unforced equation, and the term involving A as 
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the forcing or the kick. Here 5 is the usual ^-function, that is to say, the kicks occur 
instantaneously at times 0, r, 2r, 3r, . . . . 

More precisely, let $t denote the flow corresponding to the unforced equation. It 
is easy to see that for all z G §^ x M, $4(2) tends to the limit cycle 7 = = 0} as 
t —7- 00. The precise meaning of Eq. ([1]) is as follows: Let k,{6, y) = {6,y + A sin(27r6')) 
be the kick map; it represents the action of the forcing term. Then assuming we start 
with a kick at time 0, the time-r map of the flow generated by Eq. ([T]) is 

vE^^ = $^ O K, 

and the evolution of the system is deflned by iterating \E't-. We generally assume that 
r is not too small, so that during the relaxation period between kicks, the flow $i of 
the unforced equation "restores" the system to some degree. 
The parameters of interest are: 

a = amount of shear, 
A = rate of contraction to 7, 
A = amplitude of kicks, and 
r = time interval between kicks. 

Our aim in the remainder of this section is to understand - via geometric reasoning 
and numerical simulations - the meanings of these quantities, and the roles they play 
in questions of stability and chaos. Our line of reasoning follows |29j and jl4j . 

1.2 Geometry of 

A simple way to gain intuition on the geometry of is to study the \l'^-image of 7, the 
limit cycle of the unforced system. We will do so by freezing some of the parameters 
and varying others. 

Effects of varying a, A, and A 

To begin with, let us freeze A, A, and r. To flx ideas, let us take A to be relatively 
small, so that the rate of contraction is weak, and choose r large enough that e~^^ is 
a nontrivial contraction. This is when the effects of shear are seen most clearly. In 
Fig. 1, X = A = 0.1, and r = 10. Here \E't-(0,0) = (0,0) because the limit cycle 7 has 
period 1 and r is an integer multiple of this period; for non-integer r the picture is 
shifted horizontally. 

Fig. 1(b) shows four images of 7 under \I't- for increasing shear. The larger a, the 
greater the difference in velocity between two points with different y coordinates. This 
applies in particular to the highest and lowest points in /t(7) in Fig. 1(a). For a small 
enough, order in the ^-direction is preserved, i.e., for zi = (61,0) and Z2 = {ei + e,0), 
\&r(zi) will continue to have a slightly smaller 6'-coordinate than '$r{z2)- As a increases, 
some points in 7 may "overtake" others, spoiling this order. As a gets larger still, the 
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a = 0.05 



(a) The limit cycle 7 and its image ^(7) after one kick 

a = 0.25 



a = 0.5 



a = 1 





(b) ^^(7) for r = 10 



Figure 1. Effect of increasing shear. Here, A = yl = 0.1. 



total distances traveled in r units of time vary even more, and a fold develops. This 
fold can be made arbitrarily large: we can make it wrap around the cylinder as many 
times as we wish by taking a large enough. 

If we had fixed a instead, and increased A starting from A = 0.1, the resulting 
sequence of pictures would be qualitatively similar to Fig. 1(b) but in reverse order: 
The smallest A would correspond to the bottom-right image in Fig. 1(b), and the largest 
A to the top- left — provided r is scaled so that Ar remains constant. This is because 
for A small, ^(7) returns to 7 very slowly, giving the shear a great deal of time to act, 
while for larger A, ^(7) is brought back to 7 more quickly. Thus all else being equal, a 
and A, i.e., shear and damping, have opposite effects. 

The consequence of varying A while keeping the other parameters fixed is easy 
to see: the stronger the kick, the greater the difference in ^/-coordinate between the 
highest and lowest points in ^(7), and the farther apart their 6'-coordinates will be 
when flowed forward by 

What we learn from the sequence in snapshots in Fig. 1 is that A acts in concert 
with a to promote fold creation, while A works against it. 

Formulas for \E't- 

Since the unforced equation is easy to solve, one can in fact write down explicitly 
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the formulas for Let (OrjUr) = ^r(6'o7l/o)- A simple computation gives 



= eo + T+^-[yo + A sin(27r^o)] ■ (1 - e^^O (mod 1) , 

(2) 

Vr = e-^^[yo + Asm{27T9o)] . 

The reader can easily check that Eq. ([2]) is in agreement with the intuition from earlier. 

Note the appearance of the ratio jA, or rather ^(1 — e~^^), in the nonlinear term 
in the equation for 6r'. the size of this term is a measure of the tendency for a fold 
to develop in \I/t-(7). As is well known to be the case, stretch-and-fold is a standard 
mechanism for producing chaos. One can, therefore, think of the ratio 

a . shear , 

—A = ■ kick amplitude 

A contraction 

as the key to determining whether the system is chaotic, provided Ar is large enough 
that the factor 1 — e"'*'^ is not far from 1. 

Trapping region and attractor 

From the above, it is evident that much of the action takes place in a neighborhood 
around 7. Let U = {\y\ < A{e^^ - 1)"^}, so that ^r(t^) C U, and define 

F = n„>o^:^(t/) 

to be the attractor for the system Eq. ([T]). The basin of attraction of F is the entire 
cylinder §^ x M, since every orbit will eventually enter U. This usage of the word 
"attractor" implies no knowledge of dynamical indecomposability (a condition required 
by some authors). 

1.3 Lyapunov exponents 

Another measure of chaos is orbital instability, or the speed at which nearby orbits 
diverge. In this subsection, we focus on the larger of the two Lyapunov exponents of 
defined to be 

Amax(^) = lim -log = sup lim - log ■ . 

ra— i-oo n ^_^Q n— )-oo n 

Leaving technical considerations for later (see Sect. 3.1), we compute numerically Amax 
for the systems in question, sampling at various points z E U. Notice that Amax 
measures the rate of divergence of nearby orbits per kick, not per unit time. 

Each of the plots in Fig. 2 shows A^ax as a function of r for the values of a, A and 
A specified. In all six plots, we have fixed X = A = 0.1, while a varies from plot to 
plot. The first 4 values of a used in Fig. 2 are the same as those used in Fig. 1(b). 
In each plot, 10 randomly chosen initial conditions are used, the largest and smallest 
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computed values of Amax are discarded, and the largest and smallest of the remaining 
8 values are shown, the smallest as a solid black dot, and the largest, if visibly different 
than the smallest, as an open square. The plots show r G [5, 15]. We give some idea 
of the rates of contraction and sizes of the trapping regions U for these parameters: at 
r = 5, e"^^ = ^ 0.61, and U = {|y| < 0.15}; at r = 15, e'^^ = e'^-^ ^ 0.22, and 

U^{\y\< 0.03}- 

Observations from simulation results 

The discussion in Sect. 1.2 suggests that as shear is increased with other parameters 
fixed, the system is likely to get increasingly chaotic. This may lead us to expect Amax 
to increase monotonically with a. As one can see, that may be correct as an overall 
trend, but the situation is somewhat more complicated: 

In the two low-shear regimes, namely a = 0.05 and 0.25, with a few exceptions the 
computed values of A^ax are either zero or negative, with a majority of them at or very 
near zero for a — 0.05 and becoming considerably more negative at cr = 0.25. That is 
to say, Amax decreases as a increases. Notice that A^ax < means the trajectory tends 
to a sink, i.e., a stable fixed point or periodic orbit. 

Increasing shear, we see in the middle row of Fig. 2 that at first sinks dominate the 
landscape, giving way to more instances of positive Lyapunov exponents, i.e., chaotic 
behavior, as a increases. At o" = 1, the picture is very mixed, with Amax fluctuating 
wildly between positive and negative values as r varies. Notice also the nontrivial 
number of open squares, telling us that these parameters often support more than one 
type of dynamical behavior. 

In the two higher-shear regimes, a — 2 and 4, A^ax becomes more solidly positive, 
though occasional sinks are still observed. The route has been a messy one, but one 
could say that the transition to chaos is complete. 

As to the dependence on r, it appears that other things being equal, longer relax- 
ation times between kicks allow the dynamical phenomenon in effect to play out more 
completely: regardless of the sign of Amax, its magnitude increases with r in each of 
the plots. 

Finally, it is important to remember that the limit cycles used to produce the results 
in Fig. 2 are weakly attracting, making them more vulnerable to the effects of shear. 
Strongly attracting limit cycles are more robust, and larger kicks and/or shear will be 
needed to produce chaos. 

In Sects. 2 and 3 we review some rigorous theory that supports the 
numerically computed values of Amax shown. To avoid technical assumptions, we 
will focus on the model in Sect. 1.1, leaving generalizations to Sect. 4. As the reader 
will see, the mathematical ideas go considerably beyond this one example. On the 
other hand, even for this simple model, state-of-the-art understanding is incomplete. 
In the next two sections, we will vary a. A, A and r, and show that there are regions 
in the parameter space for which a clear description of the dynamics is available, and 
larger regions on which there is partial understanding. 



6 



(a) a = 0.05 



(b) a = 0.25 





fc) a = 0.5 



(d) a = 1 





(e) a = 2 





Figure 2. Lyapunov exponents Amax as functions of kick period r for increasing shear a. 
The other parameters are X = A = 0.1. For each (<J, r), we simulate 10 orbits with random 
initial conditions, iterating each for 4 x 10^ steps. We then drop the outliers and plot the 
remaining estimates, as described in the text. 
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2 Geometric Structures 



To analyze a dynamical system, it is often useful to begin by identifying its most 
prominent structures, those that are a significant part of the landscape. Even when 
they do not tell the whole story, these structures will serve as points of reference from 
which to explore the phase space. This section describes structures of this type for the 
systems defined by Eq. ([T]). 

2.1 Persistence of limit cycles at very low shear 

Proposition 2.1 ([29j) Given Ar > 0, the following hold for ^ and A sufficiently 
small: 

(a) the attractor T is a smooth, closed invariant curve near'y; 

(h) every G S"*^ x M lies in the strong stable curve W^^^{z) for some 2 G F. 

Here Wl'^{z) = {C G x M : limsup„^^ Mogrf(^;^(z), ^^^(0) < -A'r} where A' 
is a constant > |A. Proposition 12.11 follows from standard arguments in stable and 
center manifolds theory; see e.g. [ID]. The idea is simple: From Eq. (jj]), one obtains 

/ 1 + 27rf Acos(27r^)(l - e^^^ f (1 - e^^O 
\ e-^^27rA cos(27r^) e"^^ 

Since invariant cones depending on ^ and Ar clearly exist when A = 0, they will persist 
when ^ and A are small enough. 

When F is a smooth invariant curve, the dynamics on F is given by the theory 
of circle diffeomorphisms. The situation for a smooth one-parameter family of circle 
diffeomorphisms {/^} can be summarized as follows (see e.g. [7]): Let p{fu]) denote 
the rotation number of fuj. Then uj i— > p{fu)) is a devil's staircase, the flat parts 
corresponding to intervals of u on which the rotation number is rational. Moreover, 
the set of u for which p{fuj) is rational is typically open and dense, while the set of u for 
which p{fui) is irrational has positive Lebesgue measure. When p{fui) G Q, fui typically 
has a finite number of periodic sinks and sources alternating on the circle; these aside, 
every orbit converges to a periodic sink. When ^ Q, is topologically conjugate 

to an irrational rotation. 

The ideas above capture the spirit of the dynamics when shear is small enough: 
Suppose Ajnax IS computcd using an initial condition ( ^ U, and ( G W^^i^z) for 2; G F. 
Then Ajnax(C) = ^ma.x{z), and from the discussion above, the latter is either strictly 
negative or zero depending on whether p(\l/T-|r) is rational or irrational. 

Breaking of invariant curves 

Now if we fix A and A, and increase the shear a as is done in Fig. 2, the invariant 
cones - and the invariant curve itself - will break. 
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Here is how it happens in this model for integer values of r: For r G Z"*", {9,y) = 
0) is a fixed point of \1/t-, and a simple computation shows that as a increases from 

0, the larger eigenvalue of this fixed point decreases from 1 to e~^'^'^ = ^^/det{DW^, 

at which time the eigenvalues turn complex. No invariant curve can exist after that. 

Geometrically, one can think of the breaking of the invariant curve as being due to too 

much "rotation" or "twist" at this fixed point. 

Taking this observation a step further, one notes from Eq. that the rotational 

action of D'^r{0,y) is strongest at 6^ = i, where cos(27r6') = —1. This suggests that 

for fixed a, A and A, invariant curves are the most vulnerable for integer values of r, 

where this strongest rotation occurs at a fixed point. 

Interpreting Figs. 2(a) and (b) 

At 0" = 0.05, a majority of the Amax values computed are at or very slightly below 
zero. This is consistent with the existence of an invariant curve for those r. One checks 
easily that for integer values of r < 13, Amax = and the eigenvalues are complex. 

These are the first places where the invariant curve is broken as predicted. 

In the plot for a = 0.25, without pretending to account for all data points, it looks 
as though many are dropping off the A^ax = line to join the A^ax = "l-^^ liiie. The 
only holdouts for Amax = occur for smaller r where, as noted earlier, shear has not 
had enough time to act. 

2.2 Increasing shear: horseshoes and sinks 

At first, mostly sinks 

Fig. 2(c), (d) suggest that at a = 0.5, a sink with complex conjugate eigenvalues 
dominates the scene for much of the range of r considered, and the same is true at 
a = 1 for smaller values of r. 

For T G Z"*", this again is easily checked. The "twist" at ^ = | is also eminently 
visible in the last three pictures of '^ri'y) in Fig- 1(b). With a little bit of work, one can 
settle these questions rigorously, but an a priori fact that makes plausible the extension 
of this sink to non-integer values of r is that fixed point sinks with complex conjugate 
eigenvalues cannot disappear suddenly as parameters are varied: a bifurcation can 
occur only when these eigenvalues become real, i.e., ±e~2'*'^, (and a fixed point can 
vanish only when one of its eigenvalues is equal to 1). 

Finally, we remark that even though the sinks above clearly exert nontrivial in- 
fluence on the dynamics, other structures (competing sinks, invariant sets etc.) may 
be present. The many open squares in Fig. 2(c) suggest that for these parameters, 
trajectories in different regions of the phase space have distinct futures. 

Smale's horseshoes 

Horseshoes are likely present starting from a somewhere between 0.5 and 1 for r 
large enough. An example of an easily recognizable horseshoe for a = 2 and r = 10 is 
shown in Fig. 3(a). The larger a, the easier it is to give examples. 
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(b) The attractor T 

Figure 3. The attractor F and a horseshoe in it. Panel (a) illustrates the formation of 
a horseshoe: Shown are a box R (thick gray lines) and its image ^'r(-R) (thinner black 
curves). The vertical boundaries of R are the stable manifolds of the fixed point at (0, 0); 
these arc mapped into themselves by ■ Panel (b) shows a picture of the attractor T. 
The parameters are A = ^4 = 0.1, a = 2, and r = 10 for both plots. 
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Fig. 3(a) illustrates how proofs of horseshoes or uniformly hyperbolic invariant sets 
are often done: One first "spots" a horseshoe with one's eyes, namely one or more 
boxes that map across themselves in a characteristic way, and then proves that the set 
of points that remain in these boxes forever has the required splitting into expanding 
and contracting directions. When j is large, D^^^ expands strongly in the 6'-direction 
for most values of 6; contraction is guaranteed since det(D\E'^) = e^^'^ < 1. 

Proposition 2.2 ([.29J) Given \t, has a horseshoe if is sufficiently large. 

The presence of horseshoes is sometimes equated with dynamical complexity or 
chaos in the literature, and that is entirely justified insofar as one refers to the existence 
of chaotic orbits. One must not confuse the existence of these orbits with chaotic 
behavior starting from "most" or "typical" initial conditions, however: A system can 
have a horseshoe (which attracts a Lebesgue measure zero set), and have all other 
points in the phase space tending to a sink. Or, the horseshoe can be part of a 
"strange attractor", with A^ax > 0. The presence of a horseshoe alone does not tell us 
which of these scenarios will prevail. We will say more about strange attractors versus 
sinks in Sect. 3.1. Suffice it to observe here that horseshoes clearly exist for most of the 
parameters in Fig. 2(d)-(f), and A^ax is sometimes positive and sometimes negative. 

Sinks from homoclinic tangencies 

For larger shear, the attractor can be quite complicated; see Fig. 3(b). Yet in 
Figs. 2(d), (e), and even (f), sinks can also occur as noted. 

The following is purely theoretical, in the sense that we do not know exactly where 
the sinks are in these specific systems, but it is a general fact in two dimensions 
that near homoclinic tangencies of dissipative saddle fixed points ( "dissipative" means 
I det{Df)\ < 1), sinks form easily, meaning one can perturb the map and find one near 
such a tangency; see [I7j . Furthermore, tangencies persist once the stable and unstable 
manifolds of a horseshoe are shown to meet tangentially somewhere. While no results 
have been proved for this particular model, the "turns" made by unstable manifolds 
(see Fig. 3(b)) suggest the abundance of opportunities for such tangencies. 

3 A theory of strange attractors 

In this section we focus on the case of positive Lyapunov exponents, having discussed 
negative and zero values of Amax in Section 2. A combination of geometric and sta- 
tistical ideas will be used. Since these developments are more recent, we think it may 
be useful to include more background information: Sect. 3.1 discusses SRB measures 
for general chaotic systems. Sect. 3.2 surveys some recent work on a class of strange 
attractors called rank- one attractors. In certain parameter ranges, the attractors F in 
our kicked oscillator systems are of this type. In Sect. 3.3, we explain how the general 
results reviewed in Sects. 3.1 and 3.2 are applied to Eq. ([1]). 
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3.1 SRB measures 



The setting of this subsection is as follows: Let M be a Riemannian manifold or 
simply R". We consider an open set U G M with compact closure, and let / be a 
embedding of U into itself with f{U) C U. We will refer to F = r\n>of"'{U) as the 
attractor and U as its basin of attraction. Though not a formal assumption, we have 
in mind here situations that are "chaotic"; in particular, F is more complicated than 
an attracting periodic orbit. 

We will adopt the viewpoint that observable events are represented by positive 
Lebesgue measure sets, and are interested in invariant measures that reflect the proper- 
ties of Lebesgue measure, which we denote by m. For chaotic systems, the only invari- 
ant measures known to have this property are SRB measures. (The terms Lebesgue 
and Riemannian measures will be used interchangeably in this article.) 

Definition 3.1 An f -invariant Borel probability measure fi is called an SRB mea- 
sure if 

(a) Ajnax > fi-a.e., and 

(b) the conditional measures of fi on local unstable manifolds have densities with 
respect to the Riemannian measures on these manifolds. 

Recall from the Multiplicative Ergodic Theorem [T^ that Lyapunov exponents, in 
particular Amax, are defined /i-a.e., so (a) makes sense; in general, these quantities may 
vary from point to point. The meaning of (b) can be understood as follows. For an 
invariant measure to reflect the properties of m, it is simplest if it has a density with 
respect to m, but that is generally not possible for attractors: All invariant measures 
in U must live on F, and if / is volume decreasing, which is often the case near an 
attractor, then m(F) = 0. The idea of SRB measures is that if fi cannot have a density, 
then the next best thing is for it to have a density in unstable directions, the intuition 
being that the stretching of phase space in these directions leads to a smoothing of 
distributions. 

The main result on SRB measures is summarized in the next Proposition, followed 
by a sketch of its proof. The ideas in the proof will explain how SRB measures, which 
are themselves singular, are related to Lebesgue measure. Recall that for an ergodic 
measure /i, Amax is constant /i-a.e. We will denote this number by A^. 

Proposition 3.1 Let {f,^i) be an ergodic SRB measure with no zero Lyapunov expo- 
nents. Then there is a set V <ZU with m{y) > such that the following hold for every 
yeV: 

(i) Amax(l/) = A^; and 

('^'^) n Sr=o^ vif^y) ~^ I vdfJ' for every continuous observable ip : U — )■ M. 

The idea of the proof is as follows: Let 7 be a piece of local unstable manifold, and 
let be the Riemannian measure on 7. By property (b) of Definition 3.1, we may 
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assume m^-a.e. x G 7 is "typical" with respect to /i. In particular, it has properties (i) 
and (ii) in the Proposition. Let W^lx) be the stable manifold through x. Properties 
(i) and (ii) for y G W^{x) follow from the corresponding properties for x because 
d{f"'Xj f""y) ^ exponentially as n — > 00. It remains to show that the set of points 
y that are connected to /i-typical points as above has positive m-measure, and that is 
true by the absolute continuity of the stable foliation pT] . 

A little bit of history: SRB measures were invented by Sinai, Ruelle and Bowen in 
the 1970s, when they constructed for every attractor satisfying Smale's Axiom A [25] 
a special invariant measure with the properties in Definition 3.1 ( fM This 
special invariant measure has a number of other interesting properties; see e.g. [3l [31] 
for more information. At about the same time, building on Oseledec's theorem on 
Lyapunov exponents [18], Pesin [20] and Ruelle [23] extended the uniform theory of 
hyperbolic systems, also known as Axiom A theory, to an almost-everywhere theory in 
which positive and negative Lyapunov exponents replace the uniform expansion and 
contraction in Axiom A. The idea of an SRB measure was brought to this broader 
setting and studied there by mostly Ledrappier and Young; see e.g. [TT] . 

The existence problem 

While the idea and relevant properties of SRB measures were shown to make sense 
in this larger setting, existence was not guaranteed. Indeed for an attractor outside 
of the Axiom A category, no matter how chaotic it appears, there is, to this day, no 
general theory that will tell us whether or not it has an SRB measure. 

Here is where the difficulty lies: By definition, an Axiom A attractor has well- 
separated expanding and contracting directions that are invariant under the dynamics, 
so that tangent vectors in expanding directions are guaranteed some amount of growth 
with every iterate. In general, an attractor that appears chaotic to the eye must expand 
somewhere; this is how instabilities are created. But since volume is decreased, there 
must also be directions that are compressed. Without further assumptions, for most 
points X and tangent vectors v, \\Df^v\\ will sometimes grow and sometimes shrink as 
a function of n. To prove the existence of an SRB measure, one must show that on 
balance, \\Df^v\\ grows exponentially for certain coherent families of tangent vectors. 
The absence of cancellations between expansion and contraction is what sets Axiom A 
attractors apart from general chaotic attractors. 

3.2 Some recent results on rank-one attractors 

This subsection reviews some work by Wang and Young [281 SB E2] on a class of strange 
attractors. These attractors have a single direction of instability and strong contrac- 
tion in all complementary directions. Among systems without a priori separation of 

^ Sinai treated first the case of Anosov systems; his results were shortly thereafter extended to 
Axiom A attractors (which are more general) by first Ruelle and then Ruelle and Bowen. 
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expanding and contracting directions (or invariant cones), this is the only class to date 
for which progress has been made on the existence of SRB measures. 

The idea is as follows: One embeds the systems of interest in a larger collection, 
letting b denote an upper bound on their contraction in all but one of the directions 
(more precisely the second largest singular value of Dfx). One then lets 6 — )■ in what 
is called the singular limit. If this operation results in a family of well defined ID 
maps, and if some of these ID maps carry strong enough expansion, then one can 
try to conclude that for small but positive b, some of the systems have SRB measures. 
Obviously, this scheme is relevant only for attractors that have a ID character to begin 
with. For these attractors, what is exploited here is the fact that ID objects, namely 
those in the singular limit, are more tractable than the original n-dimensional maps. 

Since it is not illuminating to include all technical details in a review such as this 
one, we refer the reader to |3jy. Sect. 1, for a formal statement, giving only enough 
information here to convey the flavor of the main result: 

Let M = I X Dn-i where / is either a finite interval or the circle and -D^-i is the 
closed unit disk in R"'"^, n>2. Points in M are denoted by (x,?/) where x G / and 
y = {y^, ■ ■ ■ , y"'~^) G Dn-i, and / is sometimes identified with / x {(0, ■ ■ ■ ,0)}. Given 
F : M I, we associate two auxiliary maps: 

: M ^ M where F« = (F, 0, ■ ■ ■ , 0) , 
f : I I where / = F|7x{(o, - ,o)} • 

We need to explain one more terminology: There is a well known class of ID maps called 
Misiurewicz maps [16]. Roughly speaking, a map / is in this class if it is C^, piecewise 
monotonic with nondegenerate critical points, and satisfies the following conditions: 
(i) it is expanding away from C = {/' = 0}, and (ii) the forward orbit of every x G C 
is trapped in an expanding invariant set (bounded away from C). Maps in this class 
are known to have positive Lyapunov exponents Lebesgue-a.e. 

Theorem 1 ([31j) Let Fa : M ^ I be a 1-parameter family of maps with the 
following properties: 

(CI) there exists a* such that fa* is a Misiurewicz map; 

(C2) a ^ fa satisfies a transversality condition at a = a* and x G C{fa*); 

(C3) for every x G C{fa*), there exists j such that dyjFa*{x,0) ^ . 

Then there exists b > (depending on {Fa}) such that if Ta : M ^ M is a family of 
embeddings of M into itself with \\Ta — -f^JHc^ < b, then there is a positive measure 
set A in a-space such that for all a G A, admits an SRB measure. 

That \\Ta — F^Wqs must be sufficiently small is the rank one condition discussed 
above, and (CI) is where we require the singular limit maps to have sufficient expansion. 
Notice that (C1)-(C3) all pertain to behavior at or near a = a*. The set A will also 
be in the vicinity of this parameter. (C2) guarantees that one can bring about changes 
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effectively by tuning the parameter a, and (C3) is a nondegeneracy condition at the 
critical points. 

Remark. The existence of SRB measures is asserted for a positive measure set of 
parameters, and not for, say, an entire interval of a. This is a reflection of reality 
rather than a weakness of the result: there are parameters a arbitrarily near A for 
which Ta has sinks. In a situation such as this one where chaotic and non-chaotic 
regimes coexist in close proximity of one another, it is impossible to say for certain if 
any given map has an SRB measure. One can conclude, at best, that nearby maps 
have SRB measures "with positive probability" . 

Theorem 1 was preceded by the corresponding result for the Henon family 



The existence of SRB measures for parameters near a* = 2 and 6^1 was proved in 
^ building on results from This is the first time the existence of SRB measures 
was proved for genuinely nonuniformly hyperbolic attractors. Even though [1] is ex- 
clusively about Eq. (jl]), the techniques developed there were instrumental in the proof 
of Theorem 1. 

Returning to the setting of Theorem 1, let us call a G A a "good parameter" and 
T = Tq a "good map" . The following two properties of these maps are directly relevant 
to us. They were proved under the following additional assumption on M: 



(1) Lebesgue-a.e. z G M is contained in W^{^) where ^ is typical with respect to an 
ergodic SRB measure (in general, there may be more than one such measure). It 
follows that Ainax(-2) > for Lebesgue-a.e. z G M. 

(2) Another condition on fa* (Lyapunov exponent > log 2 and mapping every 
interval of monotonicity to all of / for some A^) implies the uniqueness of SRB 
measure. This in turn implies Amax is constant a.e. in M. 



These and a number of other results for "good maps" were proved in [32j. We 

mention one that is not used here but sheds light on the statistical properties of these 
attractors: For an SRB measure /i for which (T, yu) is mixing, the system has expo- 
nential decay of correlations for Lipschitz observables, i.e., there exists r G (0, 1) such 
that for all Lipschitz (p^tp, there exists C = C{(p,ilj) such that for all n>l, 



Ta,b ■ {x, y) ^ {1- ax^ + y, hx) . 



(4) 



det{DTa)\^V 



n-l 






15 



3.3 Application to kicked oscillators 

We now return to the model introduced in Sect. 1.1 and explain how this system can 
be fitted into the framework of the last subsection. (See [29] for details.) We fix a, A, A, 
and allow r to vary. Writing r = k + a where k = [t], the integer part of r, we let 
Tk,a = ^r- For each fixed k G Z+, we view {T^ = Tk,a,(i £ [0,1)} as the family of 
interest, and discuss if and when the conditions of Theorem 1 will hold for this family. 

Here, the singular limit maps Fa are well defined. In fact, they are the first com- 
ponents of limk^ooTk,a, i-e., 

Fa{e, y) = e + a + j-{y + A sin(27re)) , 

and the restriction of to is 

fa{e) = e + a+jAsm{27re), a E [0,1) . (5) 

Notice immediately that the range of applicability of Theorem 1 is limited to Ar 
relatively large. This is because \\Tk,a — F^Wc^ = C'(6) where b = e~^^ is required to be 
very small. For a given unforced system, where the amount of damping A is fixed, this 
means the kicks must be applied sufficiently far apart in time. 

We comment on (CI), which along with the rank one condition above are the 
core assumptions for this theorem. For our purposes let us assume fa satisfies the 
Misiurewicz condition if some iterate of fa sends its two critical points ci and C2 into 
an unstable periodic orbit or an expanding invariant Cantor set. First, such Cantor sets 
are readily available for medium size values of ^ such as ^ > 1, and unstable periodic 
orbits start to exist for somewhat smaller values of ^. Suppose for some parameter 
value a that the forward orbit of Ci is contained in an expanding invariant set K. As 
we vary a, both the orbit of ci and K will move with a. Condition (C2), assuming 
it holds, implies that for all large enough k G Z+, /^(ci) moves faster, i.e., the path 
traced out by a fa{ci{a)) cuts across K as though the latter was stationary. When 

is a Cantor set, this guarantees that /^(ci) G K for an uncountable number of 
a's. By symmetry, when that happens to /a(ci), the same is automatically true for 
/o(c2)- The larger j, the denser these Cantor sets are in S^, and the denser the set of 
parameters a that can be taken to be a* in (CI). 

The checking of (C2), (C3), and (^) are straightforward. Conditions for the unique- 
ness of SRB measures require that | be a little larger. 

To summarize, the results in the last subsection imply that for ^ > 1 (or even 
smaller), for all large enough k, there exist positive measure sets such that for 
a G Afc and r = A; + a, \E't- is a "good map" in the sense of the last subsection. In 
particular, has an SRB measure. We conclude also that Amax(-2) is well defined and 
> for Lebesgue-a.e. z E U. To ensure that Aniax(-2) = constant a.e. (so there are no 
open squares in Fig. 2) one needs to take j a little larger. It is in fact not hard to see 
that Afc ?a A/j/ for k ^ k' when both are sufficiently large, so that in this parameter 
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range, the set of r for which the properties above are enjoyed by is roughly periodic 
with period 1. 

Remarks on analytic results for chaotic systems 

Theorem 1 is a perturbative result. As is generally the case with perturbative 
proofs, the sizes of the perturbations (such as b) are hard to control. Consequently, 
applicability of Theorem 1 is limited to regimes with very strong contraction. The 
results reported in Sect. 3.2, however, are the only rigorous results available at the 
present time. Techniques for analyzing maps in parameter ranges such as those in 
Fig. 2 are lacking and currently quite far out of reach. 

The situation here is a reflection of the general state of affairs: Due to the cancel- 
lations discussed at the end of Sect. 3.1, rigorous results for the large-time behavior of 
chaotic dynamical systems tend to be challenging. 

When results such as Theorem 1 are available, however, they - and the ideas behind 
them - often shed light on situations that are technically beyond their range of appli- 
cability. Our example here is a good illustration of that: Fig. 2(d)-(f ) show that as jA 
increases, positive Lyapunov exponents become more abundant among the parameters 
tested, interspersed with occasional sinks. This is in agreement with the dynamical 
picture suggested by Theorem 1 even though with Ar G [0.5, 1.5], the contraction can 
hardly be considered strong. 

4 Generalizations 

In Sections 1, 2, and 3.3, we have focused on a concrete model. We now generalize this 
example in two different ways: 

• the unforced equation in Eq. ([1]) is replaced by an arbitrary limit cycle; 

• the specific kick in Eq. ([T]) is replaced by an arbitrary kick. 

More precisely, we consider a smooth flow $i on a finite dimensional Riemannian 
manifold (which can be M"), and let 7 be a hyperbolic limit cycle, i.e., 'y is a periodic 
orbit of period p, and for any x G 7, all eigenvalues for D^p{x) are < 1 aside from that 
in the flow direction. The basin of attraction of 7 is the set i3 := {x G M : $t(x) — )■ 7 
as t — i- 00}. We continue to consider forcing in the form of kicks, and assume for 
simplicity that the kicks are defined by a smooth embedding k : M — )■ M (as would 
be the case if, for example, k represents the result of a forcing defined by i = G{z,t) 
where the vector field G{z,t) is nonzero, or "on," for only a short time). As before, 
the kicks are applied periodically at times 0, r, 2r, 3r, . . . , and the time evolution of the 
kicked system is given by \1/t- = $t- o n. 

A new issue that arises in this generality is that one may not be able to isolate 
the phenomenon, that is to say, the kicking may cause the limit cycle to interact with 
dynamical structures nearby. Our discussion below is limited to the case where this 
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W -leaves 




Figure 4. Geometry of folding in relation to the VF**-foliation. Images of (70) ^-nd 
^' 7-2 (70) for Ti < T2, both multiples of the period of the limit cycle 7, are shown. 

does not happen, i.e., we assume there is an open set U with •y C U C B such that 
k{U) C B and $r(/t(f/)) C U, and define T = Hn^riU) to be the attractor of the 
kicked system as before. 

We further limit the scope of our discussion in the following two ways: (i) Kicks that 
are too weak will not be considered; such kicks produce invariant curves and sinks for 
the same reasons given in Sect. 2.1, and there is no need to discuss them further, (ii) 
We consider only regimes that exhibit a substantial contraction during the relaxation 
period, brought about by long enough kick intervals that permit the "shear" to act. 
As we will see, this is a more tractable situation. Rigorous results can be formulated 
- and we will indicate what is involved - but will focus primarily on ideas. Precise 
formulations of results in this generality (see [30]) are unfortunately not as illuminating 
as the phenomena behind them. 

The geometry of folding: kicks and the strong stable foliation 

As we will show, key to understanding the effects of kicks is the geometric relation 
between the kick and the strong stable foliation associated with the limit cycle of the 
unforced system. For a; G 7, we define the strong stable manifold of $t at x, denoted 
W''{x) = Wilix), to be the set W''{x) = {y e M : d($f (y), $i(a;)) ^ as t ^ cx)}; 
the distance between and in fact decreases exponentially; see e.g. [TO] . 

(This stable manifold is for the flow $j, not to be confused with that for the kicked 
map \I/t- in Prop. 2.1.) Some basic properties of these manifolds are: (i) W''^^{x) is 
a codimension one submanifold transversal to 7 and meets 7 at exactly one point, 
namely x; (ii) $((iy'^*(x)) = W^'^{^t{x)), and in particular, if the period of 7 is p, then 
^p{W^^{x)) = W^^{x); and (iii) the collection {W^^{x),x G 7} fohates the basin of 
attraction of 7. An example of a l^^'^-foliation for a limit cycle is shown in Fig. 4. 

Fig. 4 shows the image of a segment 70 of 7 under "^r- For illustration purposes, 
we assume 70 is kicked upward with its end points held fixed, and assume t = up for 
some n G where p is the period of the cycle. Since $„p leaves each ly^'^-manifold 
invariant, we may imagine that during relaxation, the flow "slides" each point of the 
curve /t(7o) back toward 7 along '^''-leaves; the larger n is, i.e., the more times it 
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laps around, the farther down the point shdes. In the situation depicted, the folding is 
quite evident. If r is not an integer multiple of p, then carries each VT^'^-manifold to 
another l^^'^-manifold. Writing t = np + a where a G [0,p), we can think of the action 
of $T- as first sliding along the iy*-manifold by an amount corresponding to $np and 
then flowing forward for time a. 

The picture in Fig. 4 gives insight into what types of kicks are likely to produce 
chaos. The following observations are intended to be informal but intuitively clear: 

(i) Kicks directed along l^^'^-leaves or in directions roughly parallel to I^'^-leaves 
are not effective in producing chaos, nor are kicks that essentially carry one W^'"^- 
leaf to another in an order-preserving fashion. For such kicks, \E't- essentially 
permutes H^'^^-leaves, and k has to overcome the contraction within individual 
leaves to create chaotic behavior. (This cannot happen in 2D.) 

(ii) The stretch-and-fold mechanism for producing chaos remains valid: the more 
"^ril) is folded, the more chaotic the system is likely to be, i.e., the intuition is 
as in Fig. 1. What is different here is that unlike our earlier example, where the 
propensity for shear-induced chaos is determined entirely by parameters in the 
unforced equation, namely a and A, we see in this more general setting that it 
matters how the kick is applied. It is the geometry of the action of k on the limit 
cycle 7 in relation to the strong stable foliation W'^'^ that determines the stability 
or chaos of the kicked system. 

(iii) The case of stronger contraction is more tractable mathematically for the follow- 
ing reason: When the contraction in is weak, as with r = ri in Fig. 4, one 
has to deal with the cumulative effects of multiple kicks, which are difficult to 
treat. When the image \E'r(7o) is pressed more strongly against 7, as in the case 
of T = T2 > Ti, cumulative effects of consecutive kicks are lessened. 

We illustrate some of the ideas above in the examples below. 

Linear shear-flow examples 

The 2D system in Eq. ([1]) : The ideas in this section can be seen as an abstraction of 
those discussed earlier. To understand that, we compute the H^^'^-leaves of the unforced 
equation in Eq. ([T]), and find them to be straight lines having slopes — ^. We concluded 
earlier that given A, the larger |, i.e., the smaller the angle between the iy**-leaves 
and 7, the more chaotic the system is likely to be. Item (ii) above corroborates this 
conclusion: Given that we kick perpendicularly to the limit cycle (as is done in Eq. ([1])), 
and points in 7 are kicked to a given height, the more "horizontal" the iy**-leaves, the 
farther the points in ^(7) will slide when we bring them back to 7. In other words, 
the J part of the ratio from earlier is encoded into the geometry of the l^^'^-foliation 
— provided that we kick perpendicularly to the cycle. 
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Generalization to n dimensions: The n-dimensional analog of Eq. ([T]) with a more 
general forcing is 

= i + ^-y, 

where 9 G S^, y G M""^, a G W^~^ is nonzero, and A is an (n — 1) x (n — 1) matrix 
all of whose eigenvalues have strictly positive real parts. For simplicity, we assume the 
kicks are perpendicular to the limit cycle {y = 0}, and to facilitate the discussion, we 
have separated the following aspects of the kick function: its amplitude is A, variation 
in 9 is H{9), and the direction of the kick is v(6') G §"~^. As a further simplification, 
let us assume v(^) = v G M"~^, i.e., it is a fixed vector. 

A computation shows that the VT^'^-manifolds of the unforced equation are given 

by 

W{9o,0) = {{9,y) -.9 = 9,- a^K-^y}, 

i.e., they are hyperplanes orthogonal to the covector {l,a^ X^^) . As noted in item 
(i) above, kick components orthogonal to (l,o""^A~^) are "dissipated" and do not have 
much effect. If if = constant, then simply permutes the H^'^^-planes and again 
no chaotic behavior will ensue. To produce horseshoes and strange attractors, a suf- 
ficient amount of variation in 9 for \E'r(7) is needed as noted in item (ii) above; that 
variation must come from H. An analysis similar to that in Eq. ([2]), Sect. 1.2, tells 
us that for large r, the amount by which the kick is magnified in the ^-direction is 
~ A H{9) cr^A^^v . We remark that the variation in H is far more important than 
its mean value, which need not be 0. 

Finally, given H ^ constant, to maximize the variation of ^E't(7) in 9 for large r, the 
discussion above suggests kicking in a direction v that maximizes cj-^A^^v. Under the 
conditions above, this direction is unique and is given by v = ±(A'^)^^o"/|(A'^)^^cr|. 
Notice that this need not be the direction with the least damping or the direction with 
maximal shear, but one that optimizes the combined effect of both. 

On analytic proofs 

When "^T- contracts strongly enough, the system falls into the rank one category as 
defined in Sect. 3.2, independently of the dimension of the phase space. This comes 
from the fact that carries all points back to 7, which is a one-dimensional object. 

If the folding (as described in item (ii) above) is significant enough for the amount 
of contraction present, then horseshoes can be shown to exist. Proving the existence of 
horseshoes with one unstable direction is generally not very difficult, and not a great 
deal of contraction is needed. 

To prove the existence of strange attractors or SRB measures, the results in Sect. 3.2 
are as applicable here as in our 2D linear example. It is proved in [30] that the periodic 
kicking of arbitrary limit cycles fits the general framework of Theorem 1, in the sense 
that as the time between kicks tends to infinity, singular limit maps are well defined. 
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They are given by /a : 7 — )■ 7, a G [0,p), where 

fa{x) := hm for all x G 7 . 

From our earlier discussion of sliding along '^-leaves, it is not hard to see that fa{x) 
is, in fact, the unique point |/ G 7 such that k{x) G W^''{y). Whether (C1)-(C3) hold 
depends on the system in question and hinges mostly on (CI), which usually holds 
when the variation is large enough. As always, these conditions need to be verified 
from example to example. 

Related Results and Outlook 

We have reviewed a set of results on the periodic kicking of limit cycles. The main 
message is that the effect of the kick can be magnified by the underlying shear in the 
unforced system to create an unexpected amount of dynamical complexity. It is an 
example of a phenomenon known as shear-induced chaos. 

A similar geometric mechanism is used to prove the existence of strange attractors in 
(a) certain examples of slow-fast systems [8] ; (b) periodic kicking of systems undergoing 
supercritical Hopf bifurcations (see [30] for details, and [I5] for results applicable to 
evolutionary PDEs); and (c) periodic forcing of near-homo clinic loops [27j. See also 
[19] . All of these results pertain to strong-contraction regimes; proofs are perturbative 
and rely on the theory of rank one attractors reviewed in Sect. 3.2. 

A welcome extension of the results reviewed here is to remove the strong-contraction 
assumption for strange attractors, but this is likely to be challenging: one has to either 
develop non-perturbative techniques or go about the problem in a less direct way. 

Random forcing is a future direction we believe to be both interesting and promis- 
ing. Numerical studies of Poisson and white-noise forcing have been carried out |14j . 
Phenomena similar to those in Sect. 1 are observed when the kick term in Eq. ([T]) is 
replaced by a term of the form 74sin(27r^) dBt where Bt is standard Brownian motion, 
i.e., the setup is a stochastic differential equation. With stochastic forcing, phase space 
geometry is messier, but A^ax depends continuously on parameters. The absence of 
wild fluctuations between positive and negative values of A^ax (corresponding respec- 
tively to strange attractors and sinks in the periodic case) gives hope to the idea that 
the analysis for stochastic forcing may be more tractable than that for periodically 
forced systems. 
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